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ABSTRACT 

The prediction of critical speeds and forced response of 
active magnetic bearing turbomachinery is of great interest due 
to the increased use of this new and promising technology. 
Calculating the system undamped critical speeds and forced 
response is important to all those who are involved in the design 
of the active magnetic bearing system. This paper is the first 
part of a two part paper which presents the theory and results of 
an investigation into the influence of sensor location on the 
undamped critical speeds and forced response of the rotor and 
bearing system. 

Part I of this paper concentrates on an extended Jeffcott 
model which was used as an approximate solution to a more 
accurate transfer matrix procedure. Theory behind a 
two-degree-of-freedom extended Jeffcott model will be presented. 
Results of the natural frequency calculation will be shown 
followed by the results of the forced response calculation. The 
system response was predicted for two types of forcing. A 
constant magnitude excitation with a wide frequency variation was 
applied at the bearings as one forcing function. The normal 
unbalance force at midspan was the second source of excitation. 
The results of this extended Jeffcott solution gives useful 
design guidance for the influence of the first and third modes of 
a symmetric rotor system. 

NOMENCLATURE 

A shaft relative motion max amplitude for 1st mode (cm) 
a mass eccentricity of imbalance (cm) 

B shaft relative motion max amplitude for 3rd mode (cm) 

C ratio of bearing damping to shaft damping (dim) 

damping of AMB (N-s/cm) 
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C 2 damping of flexible shaft (N s/cm) 

F<^ constant magnitude force applied to journal mass (N) 
i sguare root of —1, complex variable (dim) 

K ratio of bearing stiffness to shaft stiffness (dim) 

Ki stiffness of AMB (N/cm) 

K 2 stiffness of flexible shaft (N/cm) 

Kii row 1, column 1 of stiffness matrix in Jeffcott solution 

Ki 2 column 2 of stiffness matrix in Jeffcott solution 

K 21 row 2, column 1 of stiffness matrix in Jeffcott solution 

K 22 row 2, column 2 of stiffness matrix in Jeffcott solution 

L bearing span (cm) 

M ratio of bearing journal mass to rotor midspan mass (dim) 
M;il equivalent bearing journal mass (kg) 

M 2 equivalent rotor midspan mass (kg) 

R shaft absolute displacement (cm) 

X shaft maximum displacement in the X-direction (cm) 

shaft deflection at bearing journal location (cm) 

r 2 shaft deflection at midspan mass (cm) 

rg shaft deflection at AMB sensor location (cm) 

z axial distance along rotor (cm) 

Zg axial distance to AMB sensor (cm) 

a sensor relative position to midspan (dim) 

^ normalized shaft motion at sensor location (dim) 

w angular velocity of shaft (rad/s) 

0 natural frequency normalized to rigid bearing 
critical speed (dim) 
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INTRODUCTION 


The evaluation of critical speeds and forced response for 
turbomachinery with fluid-film and antifriction bearings is now 
standard design practice for many manufacturers. The standard 
transfer matrix solution technique (Myklestad, 1944; Prohl, 1945) 
is the current industry standard for evaluation of rotor response 
and undamped critical speeds. More recently interest in improved 
forced response and stability of high pressure compressors and 
pumps have forced designers to consider active magnetic bearings 
(AMBs) for either retrofit or new machinery application. The 
initial application of magnetic bearings to centrifugal 
compressor was evaluated using standard critical speed codes 
without consideration for sensor location (Hustak et al., 1987; 
Schoeneck and Hustak, 1987) . The comparisons of predicted 
response and critical speed placement to actual test and field 
results (Hustak et al., 1987; Schoeneck and Hustak, 1987; Kirk 
et al., 1988) have drawn attention to possible improvements in 
the analytical representation of the magnetic bearings. 

This paper is the first of a two part paper which presents 
an evaluation of the effect of sensor location on the predicted 
undamped critical speeds and forced response of turbomachinery. 
This paper concentrates on the solution of a 

two-degree-of-f reedom model developed by extending the original 
Jeffcott model to include bearing stiffness and damping, journal 
mass and accounting for non-colocation of bearing and sensor. 

The second paper will discuss the evaluation of a modified 
transfer matrix solution and will present results of a typical 
rotor bearing system analysis. 

The extended Jeffcott model will be considered to have 
sensors either inboard or outboard of the bearing centerline. 

The system response is calculated for two different types of 
forcing functions. The first excitation force is the usual 
unbalance located at the midspan mass. The second is a constant 
magnitude excitation applied at the journal mass while the 
excitation frequency is varied. The second type of rotor 
excitation is available in an actual active magnetic bearing and 
rotor system. 


PRINCIPLE OF ACTIVE MAGNETIC BEARING OPERATION 

The AMB is composed of two major mechanical parts, the rotor 
and the stator. Both are made of ferromagnetic laminations. The 
rotor laminations are placed on the machine shaft at the selected 
journal location. The stator laminations are slotted and include 
windings to provide the magnetic levitation and position control. 
For each degree of freedom, two electromagnets are required since 
they operate by attraction only. Figure 1 shows the stator 
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lamination construction of a radial bearing and sensor with the 
rotor lamination sleeve in the background. 


Rotor position is monitored by sensors and this signal is 
compared to a nominal reference signal with a closed loop 
controller which supplies a command signal to the . 

amplifiers. These amplifiers provide power to the 
to resist rotor movement from the nominal position. The design 
of the control loop gives the option to select the effective 
bearing stiffness and damping. The details of this design 
procedure are not the subject of this paper but the values of 
stiffness and damping must be carefully selected . 

rotor system the desired optimum dynamic response and stabii y. 

Before power is applied to the bearings, the rotor is 
supported on two auxiliary ball bearings located in close 
proximity to the AMB. The clearance between the rotor and the 
inner race of the ball bearing is selected to ppvent rotor 
contact with the AMB pole pieces or the internal seals ^he 
compressor while the rotor is at rest or during an emergency 
shutdown. When power is applied to the electronic controls, 
electromagnets levitate the rotor in the magnetic field a 
rotation by the driving source, such as a motor 

be started. The sensors and control system regulate the streng 
and direction of the magnetic fields to maintain exact rotor 
position by continually adjusting to the changing 
rotor. Should both the main and redundant features of the AMB 
fail simultaneously, the auxiliary bearing and rotor system are 
designed to permit safe deceleration. 

When the turbomachine is running the rotor shaft may take a 
dynamic mode shape such that the displacement at the sensor 
location may not be the same as the magnetic bearing centerline 
displacement. The command signal is taken from the sensor 
location but the actuator applies the force through the coil such 
that the average force acts at the bearing centerline. This 
variation in command signal and actuator location is unique o 
the active magnetic bearings and may be used to the advantage of 
the designer to help place critical speeds. The performance of 
the AMB supported machinery may be more accurately predicted it 
proper account is taken for sensor location. This requires a 
modified, iterative solution strategy for current standard 
state-of-the-art computer codes for critical speeds, forced 
response, and stability. To initially evaluate the influence of 
the sensor placement, a modified extended Jeffcott rotor model 
will be developed with an assumed deformation to study the 
sensitivity of rotor bending modes and response to sensor 
location. 
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AMB EXTENDED JEFFCOTT MODEL 


The original rotor model developed by H. H. Jeffcott consists 
of a single mass on a flexible shaft supported by rigid bearings. 
Kirk and Gunter (1972) modified this model to study the effect of 
support flexibility and damping on the synchronous response of 
the single mass flexible rotor. This paper extends the original 
Jeffcott model by assuming the existence of AMB supports. The 
extended model adds journal mass, bearing stiffness and damping 
at bearing locations, and assumes rigid bearing pedestals. The 
AMB extended Jeffcott model is shown in Figure 2. 

To develop the extended Jeffcott model the disk mass plus the 
two center quarters of the shaft mass are lumped at midspan, M 2 . 
The journal and shaft end quarter masses are lumped at bearing 
locations and modeled as Mi/2. The model is assumed to be 
symmetric; therefore, it can be simplified to a 
two-degree-of-freedom (2D0F) system as shown in Figure 3. An 
unbalance force is shown at M^, and a constant magnitude 
excitation force is shown acting on Mj^. 

The equations of motion (EOM) for the 2DOF system are written 
as follows: 


• • • • 

” 2^2 = M2w2aei‘^t _ C2(r2~rj^) - K2(r2~ri) [1] 

* • , • • • 

M^ri = F4>e^^'^ + C2(r2-ri) + K2(r2~ri) - C^rg - K^rg [2] 

In equations [1] and [2] the deflections at Mj^ and M 2 are defined 
as r^ and r 2 respectively. The deflection at the AMB sensor 
location is defined as rg. It is indicated by the EOM that the 
bearing forces are proportional to the sensor location 
deflection - not the bearing location deflection, as would be the 
case with conventional fluid-film or antifriction bearings. 

The sensor location deflection is calculated after assuming 
mode shapes of a half-period of a sine wave. These mode shapes, 
modeling the first and third modes, are shown in Figure 4. Using 
Figure 4, the equation for the sensor location deflection is 
written as follows: 

= ^1 (^2 - ^^l) sin(7ta/2) [3] 

where , 

“ = Z§/(V2) . [4] 

Equation [4] defines the value a as the ratio between the sensor 
offset and the shaft half-span. 
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After substitution of equation [3], equations [1] and [2] 
can be written in matrix form as follows: 


Ml O' 

1 

n 


'Cl o' 

‘^1 


"Ki O' 

'^1 



0 M 2 . 

1^2. 

4- 

.0 C 2 . 

.’^2. 

+ 

.0 K 2 . 

.^ 2 . 




By assuming a solution of r = the matrix equation [5] 

can be written in the following form: 


H 

1 

Ki2' 

1 

H 


Y<p 

[6] 

[k21 

K 22 . 

1^2. 





where , 

+ K2 - ) + iw ( C 2 + C2^(1-/9)) 

Ki2 = ( K^/9 - K 2 ) + iw( Ci/9 - C 2 ) 

K 21 = -Ki - iwC2 

K 22 = ( ^2 - M2W^ ) + iwC2 

p = sin Tiot/l. 


INFLUENCE OF SENSOR LOCATION ON UNDAMPED NATURAL FREQUENCY 

The influence of sensor location on the critical speeds of 
the AMB rotor system is initially investigated by calculating e 
natural frequencies of the extended Jeffcott Model. The sensor 
location is varied inboard and outboard of the bearing centerline 
by as much as 20% of half -span. 

Eliminating damping terms, and C 21 from equation [6] < and 
solving for the determinant of the resulting stiffness matrix, 
results in the following equation for the natural frequencies: 


Ki 




4 _ 


(1 - 


Ml 


(M 1 +M 2 ) 
Ml M 2 


K-,0)"^ + 


Ki K 2 


Ml M 2 


= 0 . 


[7: 


Equation [7] can be written in the following non-dimens onal 
form: 
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02 - 


[ 8 ] 


K K 

(1 (1 + K(1 - (9) )0 + - = 0 

M M 


where. 


M = M 1 /M 2 


K = K 1 /K 2 

0 = y2/(K2/M2) 


Results from Natural Frequency Analysis of AMB Extended Jeffcott 
Model 


To show how sensor location influences the first and third 
natural frequencies of various geometries of the extended 
Jeffcott model, the solution of equation [8] was graphed for 
alpha values ranging from -2.0 to 2.0. This exemplifies sensor 
separations of 20% of half-span both inboard and outboard of the 
bearing centerline. The results are shown for mass ratios, M, of 
1.0 (Figure 5), and 0.25 (Figure 6). The stiffness ratio, K, 
varies from 0.1 up to 10 in each analysis. 

The results are similar for both mass ratios. The 
sensitivity to non-colocation of bearing and sensor is increased 
in two different situations. An increase in sensitivity occurs 
as stiffness ratios increase. This is attributed to the fact 
that when the bearing stiffness increases relative to the shaft 
stiffness there is more bending energy in the rotor. This causes 
a greater difference between bearing and sensor deflection at 
higher stiffness ratios. The sensitivity also increases as 
sensor-bearing separation increases. This also results in 
greater differences between bearing and sensor deflections. 

The direction in which the criticals move depends on whether 
the sensor, at an inboard or outboard location, gives more or 
less response than the normal bearing centerline. For the first 
mode, the inboard sensors have a greater deflection than the 
bearings; therefore, the critical increases due to higher bearing 
forces. Outboard sensors detect less deflection than at the 
bearing, thus decreasing bearing forces and lowering the critical 
frequency. The opposite occurs at the third mode. The inboard 
sensors detect less deflection, while the outboard sensors detect 
more deflection than at the bearing centerline. This lowers the 
third critical frequency for inboard sensors and raises it for 
outboard sensors. 
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INFLUENCE OF SENSOR LOCATION ON FORCED RESPONSE SOLUTION 

The response is calculated for two forms of excitation 
applied independently to the AMB extended Jeffcott model 
e^ations. In order to calculate the response, the ^^^rix 
elation [6] is solved for Ri and usxng Cramer's Rule (Anton, 
1984) . The solution has the following form: 

F<#> K 22 - Ki2( ^ 2 ^^^ ) rgn 

Rl = 

Kll K22 - Ki 2 K21 

Kii( M2aw2 ) - F<p K21 

^ — — L ^ 

^11 ^22 " ^12 ^21 

The resulting response is complex, thus there exists an amplitude 
and a phase angle associated with both R^ and R 2 - The phas 
angles can be obtained from the following equations: 


0 = cos 


real (R^) 

sqrt[ {real (Ri) >2 + (imag(Ri)}2] 


02 = cos 


real (R 2 ) 

sqrt[ (real(R 2 ) + {imag(R 2 )}^l 


It must be noted at this point that the unbalance force and 
the constant excitation force are never applied at 
as this introduces additional complexities not accounted fo y 
the Jeffcott model. Physically, the constant magnitude excitation 
force is applied in one plane only, therefore the motion 
masses is in one plane only. However, 

analytical simplicity, it will be assumed that the constant 
magnitude excitation force acts in two mutually perpendicula 
plLes in a manner similar to the unbalance force. The correct 
solution for the constant magnitude excitation force ' 

then, simply the real part of the deflections R^ and R 2 shown 
above. However, since the motion is assumed to be circu ar, e 

maximum amplitudes in the X-direction, ^2' , nhase 

same as Ri and Rp , for masses and M 2 - Similarly, the ph 
angles calculated from R^ and R 2 are also valid for motion in one 
plane. Thus, regardless of whether an unbalance force or a 
Lnstant magnitude excitation force is applied, the solution 
technique and the solution itself will remain the same. 
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Results from Forced Response Analysis of AMB Extended J effnnl-t- 
Model 


results that are to be shown are of the bearing journal 
response. Similar results occur for the mid-span mass response; 
therefore, they are omitted. 

The response amplitude at the bearing location is plotted 
versus shaft frequency in Figure 7. The excitation causing this 
response is due to an imbalance at mid— span resulting from an 
eccentricity of 0.076 mm at M 2 . In this case M = 1 . 0 , K = 2 . 0 , 
and the bearing damping is set to 0.283 N-s/mm. The shaft 
damping is assumed negligible. The value of a, being varied as 
in the critical speed solution, ranges from -0.2 to 0.2. 

1*1 Figure 7 it can be seen that the first mode peak resonance 
frequencies increase from the colocation case, a = 0 , when 
inboard sensors, a > 0, are used. The peak frequencies are shown 
to decrease with outboard sensor, a < 0, use. For the third 
mode, the peak frequencies are lower with inboard sensor use, and 
for outboard sensor use. These results are very 
consistent with the results shown from the influence of sensor 
location on natural frequencies. The same reasoning can be used 
to explain both sets of results. 

The results from the case using a constant magnitude 
excitation at the bearing location shows the same tendencies as 
the unbalance case. Shown in Figure 8 are the results of the 
bearing location excitation case. 


CONCLUSIONS 

This preliminary investigation into the effect of sensor 
location on the rotor dynamic performance of AMB turbomachinery 
gives very useful results. The natural frequency and forced 
response results from the AMB extended Jeffcott model could give 
the rotor-bearing system designer greater confidence in the 
proper selection of sensor location. 

From the test run of the AMB extended Jeffcott model the 
following specific conclusions were made: 

1* For inboard sensors, as the sensor is moved away from 
the bearing, the first mode critical frequency goes higher. For 
outboard sensors, as the sensor is moved away from the bearing, 
the first mode critical frequency goes lower. 

2- For inboard sensors, as the sensor is moved away from 
the bearing the third mode critical frequency decreases. However 
for outboard sensors, as the sensor is moved away from the 


bearing the third mode critical frequency increases. 

3. The vibrational characteristics of the two mass rotor 
system did not vary when a constant excitation force was use 
instead of a rotating unbalance. 

4 The effect of the sensor position on the critical 
frequencies was considerable, when the stiffness ratio was high. 


5. Higher mass ratios led to increased sensitivity of third 
mode critical frequency to changing sensor positions. vice 
vSL? ioi bLs ritions led to m^teased deBping effects in the 
third mode, making the rotor relatively insensitive to chang g 
sensor locations. 


RECOMMENDATIONS 

The AMB extended Jeffcott model is a simple but 
very useful approximation of a much more complex rotor-bearing 
system. The following recommendations are made for 
extending the current approximate analysis: 

1. The addition of pedestal stiffness and damping should be 
included in this model. 

2 A method for accounting for sensor location should be 
included irSrtLnsfer natrlx codes that are used to calculate 
forced response. 

3 This research should be extended to a stability analysis 
in order to improve prediction capability for AMB turbomachinery. 

4. Experimental test results must be generated for 
comparison and verification of the analyses developed. 
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Figure 4 , Mode shapes for the modified JefTcott model 
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